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ABSTRACT 


To  obtain  the  expressions  for  mean  and  variance  of  reliability  estimate,  R = P(X  > Y)  , analytically,  is  generally 
difficult.  Here,  we  find  approximate  expressions  for  mean  and  variance  of  estimated  system  reliability  in  interference 


theory  when  stress  and  strength  follow  some  particular  distributions.  We  have  evaluated  approximates  of  mean  and 
variance  of  estimated  reliability  when  stress-strength  both  follows  either  exponential  or  normal  distribution.  For  validity  of 


approximation  method,  we  have  used  Monte-Carlo  simulation.  Also,  Normal  probability  plots  of  estimated  reliability 
samples  are  drawn  for  different  values  of  the  parameters  of  the  distribution.  From  Monte  Carlo  simulation  (MCS),  it  is 
observed  that  approximation  for  mean  and  variance  of  estimated  reliability  is  up  to  the  mark. 

KEYWORDS:  Reliability,  Interference  Theory,  Monte-Carlo  Simulation  (MCS) 

1.  INTRODUCTION 

In  Interference  theory  of  reliability  R=P(X>Y),  the  reliability  of  a system  and  its  other  reliability  characteristics 
can  be  expressed  as  some  functions  of  the  parameters  of  the  distributions  of  strength  (X)  and  stress  (Y)  associated  with  the 
functioning  of  the  system.  We  estimate  these  parameters  and  substitute  these  values  in  the  expressions  for  reliability  and 
other  characteristics  to  get  their  estimates.  If  the  estimates  of  parameters  used  here  are  maximum  likelihood  estimators  then 
from  the  invariance  property  of  MLE’s,  the  corresponding  estimators  of  reliability  are  also  MLE’s.  There  exists  extensive 
literature  for  estimation  of  reliability  analytically  for  single  component  systems.  But  the  reliability  expressions  for 
multi -component  systems  are  not  simple  enough  to  facilitate  analytical  estimation  of  reliability  and  its  other  characteristics. 
Also,  due  to  lack  of  stress-strength  data  one  way  out  is  simulation.  For  example,  (Manders  et  ah,  1982,  Aldrisi,  1987, 
Stumpf  and  Schwartz  1993,  Zhang  et  ah,  2010)  have  simulated  stress-strength  and  estimated  reliability.  (Paul  and 
Borhanuddin,  1997,  Rezaei  et  ah,  2010)  estimated  reliability  of  stress-strength  model,  using  Monte-Carlo  simulation 
(MCS).  (Ahmad  et  ah,  1997),  obtain  Bayes  estimates  of  P(Y<X)  using  MCS.  (Uddin  et  ah,  1993)  estimated  reliability  for 
multicomponent  system  using  MCS.  (Patowary  et  ah,  2012)  estimated  reliability  of  n-standby  system  using  Monte-Carlo 
simulation.  Similarly,  it  is  difficult  to  obtain  the  distribution  of  a reliability  estimator  or  even  its  exact  mean  and  variance, 
analytically.  In  this  paper,  an  attempt  has  been  made  to  find  approximate  expressions  for  mean  and  variance  of  estimated 
reliability  on  the  basis  of  (Lyold  and  Lipow,  1962)  and  check  the  validity  of  the  approximation  by  Monte-Carlo  simulation. 

In  Section  2,  we  have  given  the  approximate  expressions  for  mean  and  variance  of  single  parameter  and  two 
parameters  distribution  cases.  In  Section  3,  approximation  expressions  of  mean  and  variance  are  obtained  when  stress- 
strength  (S-S)  both  follows  exponential  or  normal  distribution.  Also,  Monte-Carlo  simulation  (MCS)  is  extensively 
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performed  for  verification  purpose  of  expressions. 

2.  METHOD  TO  OBTAIN  MEAN  AND  VARIANCE  OF  ESTIMATED  RELIABILITY 

In  this  paper,  existing  method  of  (Lloyd  and  Lipow,  1962)  is  used  to  estimate  reliability  in  interference  models 
when  stress-strength  follows  particular  distribution.  However,  as  it  is  not  easy  to  get  real  life  data,  we  have  used  MCS  for 
estimation  of  mean  and  variance  of  estimated  reliability. 

The  reliability  can  be  expressed  as  a function  of  the  parameters  0 = 0[,  0->,  ...  of  the  distributions  of  strength 

X and  stress  Y.  According  to  the  above  method,  if  R ( 0 ) is  the  MLE  of  R ( 0 ),  then 

(i)  If  0 = 0 i.e  only  one  parameter  is  involved,  then  the  approximate  mean  of  R(0)  is  given  by 


ZV  /V 

rn 

E 

R(0) 

= R (0)  + oj 

Im  J 

when  E ( 0 ) ~ 0 , 


Where  O 


vM, 


-AO  as  M ->  OO  . 


(2.1) 


M is  the  sample  size  of  estimated  reliability  samples.  Obviously,  R(0)  is  asymptotically  unbiased  for  R(  0 ). 


The  approximate  variance  of  R(0)  is  given  by 


r n 

zv  zv 

3R(0) 

A 

f 1 ^ 

Var 

R(0) 

00 

Var  ( 0 ) + 0 

9=9 

where  O 


' 1 ^ 

vM3/2J 


->o  as  M -A  oo  . 


(ii)  If  0=  (X,|T),  i.e.,  two  parameters  case,  then  the  approximate  mean  and  variance  of  R ( X , p ) are  given, 
respectively,  by  Eq.(2.3)  and  Eq.(2.4) 


R(X,  p)  =R(X,|t)  + 0 


V My 


(2.3) 


when  E( X)  ~ X and  E( p)  ~ (I . 

So,  R (X,  p)  is  asymptotically  unbiased  for  R(  X,  (J.) . 


And  Var 


R(X,  p) 


3R(A,,P) 

dX 


~l2 


Var(X)  + 


-I  A, = A,, |i=|»i 


3R(AP) 
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Var(p) 


l=>.,p.=p. 
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Let  DM  be  the  difference  between  E 


DM  = 


Cov(A.,(l)  + 0 

A.=A.,p.=pt 

R(X,  (I)  and  R(  Al,  |I)  i.e. 


1 


m 3/2 ; 


(2.4) 


R(X,  A)J-R(l,|i) 

Similarly,  let  DV  be  the  difference  of  V ar 


(2.5) 


R(Ai,  A) 


and  r.h.s  of  Eq.(2.4). 


If  the  values  of  these  differences  are  negligible,  then  approximate  expressions  are  considered  as  satisfactory. 

3.  STRESS-STRENGTH  FOLLOWS  PARTICULAR  DISTRIBUTION 

We  have  seen  that  in  interference  models  system  reliability  is  a function  of  stress-strength  parameters.  Let  f(x)  be 
the  p.d.f.  of  strength  (X)  of  the  system  and  g(y)  be  that  of  the  stress  (Y)  on  the  system.  Here,  we  have  considered  two 
cases. 

Case  I:  When  both  f(x)  and  g(y)  follows  exponential  distribution  (one  parameter  case) 

Case  II:  When  both  f(x)  and  g(y)  follows  normal  distribution  (two  parameter  case) 

3.1  Stress-Strength  Exponentially  Distributed 

Let  f(x)  and  g(y)  be  exponential  with  means  A.  and  |i,  respectively.  Then  the  system  reliability  is  given  by 
(Kapur  and  Lamberson,  1977) 


A,  . . 

R = = R ( A,, (I),  say.  A,  > 0,  (J,  > 0 

A.  + [X 


(3.1.1) 


Suppose  M units  are  put  on  test.  Let  xi,  x2,  ...,  xM  be  the  strengths  of  the  M units  and  let  yi,  y2,  ...,  yM  be  the 
stresses  working  on  them.  Now  we  know  that  for  exponential  distribution  sample  mean  is  an  MLE  of  population  mean,  and 


is  unbiased,  consistent,  sufficient  and  efficient.  Hence,  X 
the  same  properties, 

i.e.  A.  = X and  A = y 
An  estimate  of  R is  given  by 


f 1 V 1 

\ 

II 

M 

and  y 

l M J 

^ M J 

are  MLE’s  of  Aland  (J,  with 


(3.1.2) 


R = 


x + y 


= R(x,  y)  , say 


(3.1.3) 


Since,  R is  a one-valued  function  of  X and  y,  hence  from  the  properties  of  MLE’s,  R is  an  MLE  of  R.  So 
from  Eq.(2.3) 
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( ] \ 

E R(x,  y)  = R(X,  |J.)  + 0 — , 

J M 


(3.1.4) 


Thus  R(x,  y)  is  asymptotically  unbiased  for  R(X,  (i) . 


Further,  since  X and  Y are  independent  hence  X and  y are  also  independent  and  so 


Cov(x,  y)  = 0 


Then  from  Eq.(2.4)  and  Eq.(3.1.5)  we  have 


(3.1.5) 


’R(x,  y)"  = 


3R(x,y) 

3x 


Var(x)  + 


x=X,y=p 


0R(x,y) 


Var(y) 


x=>.,y=p 


(3.1.6) 


Now,  from  Eq.(3.1.3) 


0R(x,y)  = Y 

3x  (x  + y)2 


9 R(x,  y) 


(x  + y)2 


(3.1.7) 


We  also  know  that  for  exponential  distributions,  considered  here 


Var  (X)  = X2  and  Var(Y)  = (T 


Hence,  Va,(X,  . YE™. - + and  Van,  Y , . Y!£<X> . HI 
MM  MM 

Substituting  from  Eq.(3.1.7)  and  Eq.(3.1.8)  in  Eq.(3.1.6),  after  some  simplifications,  we  get 


(3.1.8) 


Var  R(x,  y)  = 


2 X2[L2  { 1 

M(X  + \i)4  {.M31 


(3.1.9) 


which  ->0  as  M->  °° . Hence,  R(x,  y)  is  consistent  for  R(X,  |4) . Thus  R(x,  y)  is  asymptotically  unbiased, 
consistent  and  sufficient  for  R(X,  |J,)  , since  X,  y are  sufficient  for  X and  (4,  respectively. 

From  Eq.(3.1.4)  and  Eq.(3.1.9)  we  get 


DM  = E 


R(x,y)  -■ 


(3.1.10) 


DV=  Var 


R(x,  y)  ■ 


2 X2\i2 

M(A,  + n) 


(3.1.11) 
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We  use  MCS  to  establish  the  validity  of  the  approximate  expressions  Eq.(3.1.4)  and  Eq.(3.1.9).  As  earlier,  we  take 
independent  samples  of  sizes  M from  exp(X)  and  cxpffl) , populations  for  a given  X and  (J,.  From  these  we  obtain 

estimates  of  X and  p as  X and  y and  substituting  these  values  in  Eq.(3.1.4)  we  get  an  estimate  of  R . The  whole  process 

is  repeated  k times  and  accordingly  we  get  k values  of  R . We  obtain  mean  and  variance  of  R viz.  E[R(x,  y)]  and  Var 

A /X 

[R(x,  y)] . The  values  are  tabulated  in  Table  3.1.1.  We  have  also  drawn  normal  probability  plot  (NPP)  for  R in  each  case 
and  found  that  the  fitting  of  normal  distribution  is  quite  good.  It  conform  finding  of  (Lloyd  and  Lipow,  1962).  For 
illustration  purpose,  we  have  given  a NPP  graph  for  X=  2,  p=  2,  M = 500  and  k = 100  in  Fig.3.2.1.  DM  and  DV  are 
obtained  from  Eq.(3.1.10)  and  Eq.(3.1.11)  for  given  X,  (J,  and  M and  given  in  the  same  Table  3.1.1. 


Figure  3.1.1:  (Exponential  S-S) 


Note:  In  Table  3.1.1,  V = 


1 2xy- 


Table  3.1.1:  Exponential  Stress-Strength 


M 

K 

X 

R 

Mean  of  R 

SD  of  R 

V 

Z 

DM 

DV 

500 

50 

2 

2 

0.5000 

0.5023 

0.0139 

0.0500 

1.3014 

0.0006 

0.0023 

500 

200 

2 

2 

0.5000 

0.4990 

0.0152 

0.0250 

0.9045 

0.0007 

0.0004 

500 

100 

2 

2 

0.5000 

0.4992 

0.0155 

0.0354 

0.5232 

0.0012 

0.0011 

500 

100 

3 

2 

0.6000 

0.6017 

0.0156 

0.0339 

0.0827 

0.0011 

0.0009 

500 

100 

2 

3 

0.4000 

0.3993 

0.0150 

0.0339 

0.4641 

0.0016 

0.0009 

100 

100 

2 

2 

0.5000 

0.4962 

0.0390 

0.0354 

0.9792 

0.0045 

0.0014 

100 

100 

2 

3 

0.4000 

0.3995 

0.0328 

0.0339 

0.1573 

0.0041 

0.0014 

100 

100 

3 

2 

0.6000 

0.5996 

0.0290 

0.0339 

0.1372 

0.0001 

0.0012 

50 

100 

2 

2 

0.5000 

0.5029 

0.0551 

0.0354 

1.6772 

0.0019 

0.0001 

50 

100 

2 

3 

0.4000 

0.4097 

0.0478 

0.0339 

2.0269 

0.0006 

0.0001 

50 

100 

3 

2 

0.6000 

0.5050 

0.0562 

0.0339 

0.8949 

0.0078 

0.0004 

As  noted  earlier,  for  achieving  a better  estimates  of  R,  we  have  taken  M = 500.  But  to  see  the  accuracy  of 
Eq.(3.1.4)  and  Eq.(3.1.9),  we  have  taken  M = 50  and  100  also  and  found  that  all  z-values  are  insignificant  barring  M = 50 
when  k = 100,  X=  2 and  (1=  2.  The  values  of  DM  and  DV  are  quite  small  for  different  values  of  M,  so  we  may  assume 
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that  approximations  given  by  Eq.(3.1.4)  and  Eq.(3.1.9)  are  up  to  the  mark. 

3.2  Stress  -Strength  Are  Normal  Variates 

Let  us  suppose  that  X~N(  A.,  GT  ) and  Y~N(  JI,  G2).  Then  we  have  seen  that  the  reliability  of  the  system  is  given 
by  (Kapur  and  Lamberson,  1977) 


R = 


O 


V(032  + G2) 


- °o  < X,  |I  < oo  ; > 0,  G > 0, 


(3.2.1) 


where  <£(*) 

is  c.d.f.  of  standardized  normal  variate. 

Without  loss  of  generality,  we  may  assume  that  one  of  the  variable  follows  N(0,  1).  Or  often  the  mean  and 
variance  of  strength  are  known  because,  eg.,  production  is  more  or  less  under  control.  In  such  case  also  we  may  assume 
that  X ~ N(0,  1).  Then  Eq.(3.2.1)  gives. 


R = 0 


= R(|T,  G")  , say 


(3.2.2) 


Let  us  suppose  that  M components  are  put  on  test  whose  strengths  are  known  and  let  yi,  y2,  ...,  yM  be  the  M 

— 2 2 

stresses  on  them.  We  know  that  y and  s are  independent  MLE’s  of  (I  ando  , respectively,  which  are  unbiased, 
consistent  and  sufficient,  where 


s 


2 


g(yj-y)2 

i=i  M-l 


l M 

and  y =77  £ yi 
M i=i 


(3.2.3) 


Then  from  the  properties  of  MLE’s,  R , given  below,  is  the  MLE  of  R 


R = 1-0 


V(  1 + s2) 

Hence,  from  Eq.(2.3) 


R(y,  S~) , say 


R(y,  s2)  =R(n,  oz)  + o 


f i T 


vMy 


Now,  since  y and  s2  are  independent,  hence 


Cov(y,  S2)  = 0.  (3.2.6) 

Hence  from  Eq.(2.4)  and  Eq.(3.2.6)  we  have 


(3.2.4) 


(3.2.5) 
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3R(y,s2) 

Var 

R(y,s2)  = 

dy 

— 2 0 

y=p,s  =<r 


Var(y)  + 


8R(y,s  ) 
ds2 


Var(s2) 


— 2 2 
y =(x,s  =o 


+ o 


f 1 > 

i3/2 


vM  y 

Now  from  Eq.(3.2.2) 


3R(y,s2) 

1 

dy 

Vd  + s2) 

$ 


Vd  + s2) 


3R(y,s2) 

1 

y 

<P 

y 

ds2 

2 

[_(l  + s2)3/2  J 

_Va+s2)_ 

And 


Further,  we  know  that 


a2  2 2a4 

Var(  y ) = — and  Var(s')  = 

M M 


Substituting  from  Eq. (3. 2. 8) , Eq.(3.2.9)  and  Eq.(3.2.7)  in  Eq.(3.2.10)  we  get 

1 a2  2 M2  2a4  J 1 A 

O + 2 3 — - <£  + O 


Var 


R(y,  s2) 


Where  $ = $ 


1 + a2  M 4(l  + o-)J  M 

“M- 


m3/2; 


V(i+°“) 


Or  Var 


R(y,  s2)' 


Q2 

M 


_2  2 4 

O LI  O 

- + ^ 


l + o2  2(1 + o2)3 


(3.2.7) 


(3.2.8) 


(3.2.9) 


(3.2.10) 


(3.2.11) 


(3.2.12) 


Thus,  we  see  from  Eq.(3.2.5)  that  R is  asymptotically  unbiased  for  R and  from  Eq.(3.2.11)  it  is  consistent  also. 

— 2 1 
Since,  y , s are  sufficient  for  (J.  and  (7" , so  it  is  sufficient  also. 


Now,  from  Eq.(3.2.5)  and  Eq.(3.2.1 1), 
E 


DM= 


And  DV  : 


R(y,s2)  -R(fi, a2) 


Var[R(y,s2)]-  ' 


a2  2 


2a  2 


(3.2.13) 


(3.2.14) 


l + o2  M 4(1  + a2)3  M 

Here,  also  we  use  MCS  for  the  validation  of  the  approximations  Eq.  (3.2.5)  and  Eq.(3.2.11).  For  this  as  in  the 
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Sec. 3.1,  we  take  samples  of  different  sizes  M from  N(|l,  0 ) population  for  particular  values  of  (J,  and  0" ; its  mean  and 


variance  are  y and  s2,  respectively.  Substituting  these  values  in  Eq.(3.2.4),  we  get  an  estimate  R of  R.  The  whole  process 

is  repeated  k times  giving  k values  ofR  . Its  mean  and  variance  gives  E[R(y,  S~)]  and  Var[R(y,  S^)],  respectively. 
DM  and  DV  are  obtained  from  Eq.  (3.2.13)  and  Eq.(3.2.14).  The  values  are  tabulated  in  Table  3.2.1.  We  have  also  drawn 
NPP  graphs  for  R in  each  case  of  Table  3.2.1.  NPP  graphs  suggest  that  the  data  sets  reasonably  follow  the  normal 
distribution.  For  illustration  purpose  we  have  given  only  one  NPP  graph  of  R when  M = 500,  k = 200,  [1=1  and  0=1 
in  Figure  3.2.1. 

Here,  also  we  have  taken  M = 500,  but  for  checking  the  accuracy  of  Eq.(3.2.5)  and  Eq.(3.2.11),  we  have  taken  M 
= 50  and  100  also.  From  Table  3.2.1,  we  note  that  for  k = 100,  z-value  is  significant.  So  for  achieving  a better  estimate  of  R 
we  have  taken  k = 200  throughout  and  found  that  all  z-values  are  insignificant.  For  example,  z = 0.2522  when  m = 0 and 
0=1.  We  observe  that  the  values  of  DM  and  DV  are  quite  small  for  different  values  of  M,  so  the  approximations  given  by 
Eq.  (3.2.5)  and  Eq.  (3.2.11)  are  reasonable. 


Table  3.2.1:  Normal  Stress-Strength 


M 

K 

ft 

o 

R 

Mean  of  R 

SD  of  R 

Z for  R 

DM 

DV 

500 

100 

0 

1 

0.5000 

0.4975 

0.0115 

2.2126 

0.0002 

0.0001 

500 

200 

0 

1 

0.5002 

0.5002 

0.0122 

0.2522 

0.0002 

0.0001 

500 

200 

1 

1 

0.2398 

0.2400 

0.0116 

0.2822 

0.0002 

0.0001 

500 

200 

2 

1 

0.0787 

0.0787 

0.0056 

0.1737 

0.0000 

0.0000 

500 

200 

-1 

1 

0.7603 

0.7608 

0.0113 

0.7002 

0.0007 

0.0006 

500 

200 

-2 

1 

0.9214 

0.9211 

0.0057 

0.6402 

0.0001 

0.0011 

100 

200 

0 

1 

0.5000 

0.4986 

0.0277 

0.7397 

0.0016 

0.0010 

100 

200 

1 

1 

0.2398 

0.2403 

0.0252 
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Figure  3.2.1:  (Normal  S-S) 
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NAAS  Rating  3.45 
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4.  CONCLUSIONS 

In  this  paper,  approximate  expressions  for  mean  and  variance  of  estimated  reliability  in  interference  models  are 
evaluated  using  (Lloyd  and  Lipow,  1962)  method.  Here,  we  have  considered  two  cases-when  both  stress  and  strength 
follow  exponential  distribution  and  when  both  stress  and  strength  follow  normal  distribution.  Validity  of  approximations  is 
checked  by  Monte-Carlo  simulation.  We  have  found  that  expressions  are  acceptable,  particularly  when  the  sample  size  is 
reasonably  large..  The  approximate  expressions  can  be  extended  when  stress  strength  follow  more  than  two  parameters 
distributions,  also. 
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